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Abstract 

A model of simplicial quantum gravity in three dimensions is investigated numer- 
ically based on the technique of the dynamical triangulation (DT) . We are concerned 
with the surfaces appearing on boundaries (i.e., sections) of three-dimensional DT 
manifold with 5 3 topology. A new scaling behavior of genus distributions of bound- 
ary surfaces is found. Furthermore, these surfaces are compared with the random 
surfaces generated by the two-dimensional DT method which are well known as a 
correct discretized method of the two-dimensional quantum gravity. 
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1 Introduction 



There has been, over the last few years, remarkable progress in the quantum theory of 
two-dimensional gravity. Two distinct analytic approaches for quantizing two-dimensional 
gravity have been established. These are recognized as a discretized JlJ and a continu- 
ous p| theory. The discretized approach, implemented by the matrix model technique, 
exhibits behavior found in the continuous approach, given by Liouville field theory, in a 
continuum limit. It thus seems to exist strong evidence for the equivalence of the two 
theories in two dimensions. Numerical methods based on the matrix model, such as that 
of dynamical triangulation || have drawn much attention as alternative approaches to 
studying non perturbative effects, being also capable of handling those cases where an- 
alytical theories cannot yet produce meaningful results. In the dynamical triangulation 
method, calculations of the partition function are performed by replacing the path inte- 
gral over the metric to a sum over possible triangulations. Several close studies [|J on the 
structure of two-dimensional DT surfaces in quantum gravity were made, and it revealed 
fractal structures of these surfaces. 

However, we feel that the relationship between three-dimensional (also four-dimensional) 
quantum gravity and its discretized model is not yet understood. Over the past few years 
a considerable number of numerical studies have been made on three-dimensional simpli- 
cial quantum gravity |J. It is generally agreed that the phase transition of three- and 
four- dimensional simplicial quantum gravity are both first order |||7|. The problems 
are still in controversy. On the other hand, recent numerical results obtained by the dy- 
namical triangulation for three and four dimensions suggest the existences of the scaling 
behavior near to the critical point Hl^llOl- It may be important to note the asymmetric 
behavior of the order parameter in usual phase diagram. When the coupling strength, n , 
closes to the critical point from the strong coupling side (k < Kq) which corresponds to 
the crumple phase, the transition is smooth. On the other hand, when kq closes to the 
critical point from the weak coupling side (k > Kq) which corresponds to the branched 
polymer phase, the transition is very rapid. In ref. |J it is reported that near the critical 
point belonging to the strong coupling phase the scaling behavior of the mother universe 
exists and also that no mother universe exists in the weak coupling phase (i.e., branched 
polymer phase). It seems reasonable to suppose that the model makes sense as long as 
we close to the critical point from the strong coupling side. We actually observe a dou- 
ble peak histogram structure which is a signal of the first order phase transition for an 
appropriate large lattice size. Therefore, we carefully chose one of peaks belonging to the 



strong coupling phase as an ensemble for our simulations [ 1 



This paper is intended as two investigations, firstly, the searching a new scaling prop- 
erty for the boundary surfaces in three-dimensional DT manifolds near the critical point 
(kq) and, secondary, the relation between these boundary surfaces and two dimensional 
random surfaces implemented by the matrix model. Below we will consider only three- 
dimensional DT manifolds with S 3 topology. We can represent a typical configuration in 
Fig.[I]. M denotes a three-dimensional DT manifold with S 3 topology and Si, S2 and S3 
denote the boundaries which are closed and orientable triangulated surfaces at a distance 
d from an origin (x). The geodesic distance can be recognized as a usual time parameter. 

Our model of three-dimensional Euclidean quantum gravity with boundaries may be 
recognized as a model of a quantum nucleation of the universe in (2 + l)-dimensional 
gravity. The nucleation of the universe by a quantum tunneling may be described by 




Figure 1: Schematic picture of boundary surfaces (S 1; S 2 and S 3 ) at distance d and a 
surface (E 2 ) & t distance d! in 3D Euclidean space M with S 3 topology. The mother 
boundary surface is defined as a surface (E 2 ) with the largest tip volume (V 2 ), and the 
other surfaces are defined as the baby boundary. 



going out of the Euclidean signature region to the Lorentzian signature region in a sense 



of the semiclassical approximation |TT|. The nucleation of the universe can be regarded 
as a topology-changing process in the sense that the universe takes a transition from the 
initial state with no boundary to the final state with nontrivial topology. We can treat 
the model in a full quantum way which means that we sum up all of the fluctuations of 
the three-dimensional metric g^ and also put no restriction on the boundary surfaces E. 
Of course, it is possible to introduce the extrinsic curvature on the boundary surface as a 
physical restriction. The process of a quantum tunneling requests that all the components 
of the extrinsic curvature vanish (i.e., totally geodesic). 

This paper is organized as follows. In Sec. 2 we briefly review the model of the three- 
dimensional dynamical triangulation. In Sec. 3 all of our numerical results, especially genus 
distributions and coordination number distributions for the mother boundary surfaces, are 
shown. We summarize our results and discuss some future problems in the final section. 



2 Model 

It is not yet known how to give a constructive definition of three-dimensional quantum 
gravity. We start with the Euclidean Einstein-Hilbert action, 

S EH = J d 3 x^ (a - ifl) , (1) 

where A is the cosmological constant, and G is Newton's constant of gravity. We use the 
lattice action of the three-dimensional model with the S 3 topology, corresponding to the 
above action, as follows: 



S(k ,k 3 ) = -k N + K3N3 



= -^iV +(A'-^(27r-6cos- 1 (i)))iV3, (2) 

where Aj denotes the total number of i-simplexes, and A' = cA; c is the unit volume, and 
cos _1 (|) is the angle between two tetrahedra. The coupling kq is proportional to the in- 
verse of bare Newton's constant, and the coupling n 3 corresponds to a lattice cosmological 
constant. 

For the dynamical triangulation model of three-dimensional quantum gravity, we con- 
sider a partition function of the form 

Z( Kq ,k 3 ) = ]T e" 5 ^). (3) 

T(S 3 ) 

We sum over all simplicial triangulations, T(S 3 ), on a three-dimensional sphere. In prac- 
tice, we must add a small correction term AS to the lattice action in order to suppress 
volume fluctuations. The correction term is denoted by 

AS = 5(N 3 - N^ taTget) ) 2 , (4) 

where iV^ ax S e ) - IS ^ e target value of three-simplexes, and we use 5 = 0.0005 in all our 
run. A fine tuning of k,q is one of the elaborate works in this model because the critical 
coupling Kq depends on the size of the system. A Ko^ U1V discussed in ref. |L2| is used as 



the critical value Kq(N^) in our simulation. In ref. |12|] the parameter kq is tuned to Kq^ U1V 
using the multicanonical Monte Carlo method such that the heights of double peaks of 
Nq become the same. 



3 Measurements 

We now define the intrinsic geometry using the concept of a geodesic distance as a mini- 
mum length (i.e., minimum step) in the dual lattice between two tetrahedrons in a three- 
dimensional DT space. Suppose a tree-dimensional ball (3-ball) which is covered within 
d steps from a reference 3-simplex in the three-dimensional manifold with S 3 topology. 
Naively, the 3-ball has a boundary with spherical topology (S 2 ). However, because of 
the branching of the DT space, the boundary is not always simply-connected, and there 
usually appear many boundaries which consist of closed and orientable two-dimensional 
surfaces with any topology and nontrivial structures such as links or knotsQ. These two- 
dimensional boundary surfaces are equivalent to two-dimensional randomly triangulated 
surfaces with corresponding genuses. In order to discuss the scaling properties of these 
surfaces, the genus distributions of these surfaces are measured. We should notice that 
there appear many boundaries with distance d (see Fig.|lJ) . The boundaries are divided 
into two classes: one is a baby boundary and the other is a mother one. The boundaries 
with small sizes (~ 0(1)) are called by a baby one. The baby boundaries are originated 
from the small fluctuations of the three-dimensional Euclidean spaces. We thus think 
that these surfaces are non-universal objects and then suffer from large finite size effects. 
In this section, we shall concentrate only on the mother boundary. 

Here, we give a precise definition of "baby" and "mother" universes in Fig.|I]. The 
mother universe is defined as a boundary surface (E 2 ) with the largest volume (V2), 



(a) 



(b) 



Figure 2: Measurement of a genus of boundary surface, (a) denotes an example of a 
irregular configuration. The boundary surfaces in general may have many singular vertices 
and links, (b) denotes a regular configuration. After we change (a) into (b), we measure 
the genus of (b). 

and the other surfaces are defined as baby universes. The problem is whether we can 
distinguish these two kinds of universes (baby and mother boundary surfaces) in our 
simulations. The answer is yes. In Fig.|3| we plot the tip volumes (correspond to V\, Vj> 
and V3 in Fig.l) distributions near to the critical point (k = 4.195 and k 3 = 2.222) with 
appropriate geodesic distances when jV 3 = 32. This figure tells us that there is one mother 
boundary whose tip volume is prominent and the others belong to the class of the baby 
universe. At any rate, we can easily distinguish the mother boundary surface from the 
others by the measurement of the tip volume. 

In ref. [§] it is reported that the surface-area-distributions (SAD) of the mother uni- 
verse show the scaling behavior near to the critical point but the baby universe does not 
show a scaling behaviour at all. The important point to note is that in ref. || the mea- 
surements in the critical region were done in the strong coupling phase. Therefore, we 
focus on the mother universe near to the critical point taking the limit from the strong 
coupling phase following subsections 3.1 and 3.3. On the other hand, it is known that in 
the weak coupling limit the three-dimensional DT manifold becomes a branched polymer, 
and its boundary surfaces show no scaling property ||. Therefore, we cannot consider 
the genus distributions! of boundary surfaces in the weak coupling phase. 

3.1 Genus distributions of boundary surfaces near to the crit- 
ical point 

An interesting observable for three-dimensional DT manifolds is a genus distribution of 
boundary surfaces with various geodesic distances. The genus of the mother boundary 
surfaces is, in fact, naively expected to scale because the genus is a dimensionless quantity. 
The boundary surfaces in general have very complicated and irregular structures. We 
illustrate an irregular surface as a simple example: (a) in Fig|| In order to calculate the 
genus of these boundary surfaces we must dispose these irregular configurations. We thus 
detach irregular sub-simp lexes (vertices or links) like in Fig.0 ((a) — > (b)) i 

*In the strict sense links or knots are constructed by loops. In our case the loop is a fat loop. 

^We obtain only spherical boundary surfaces at any geodesic distances in this phase. 

■''This deformation ((a) — > (&)) is recognized as the infinitesimal deformation d — > d — e. Since e can 
be arbitrary small, we obtain the same distributions of the boundary genus as at distance d in the limit 
£ -> 0. 




Figure 3: Typical distributions of the tip volume distributions for all the boundaries with 
N 3 = 32JT and d — 10, 11 and 12. Discrepancies shown in a small part of the tip volume 
(10 ~ 100) are due to the finite size effects. 
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Figure 4: Genus distributions of mother boundary at the critical point for various dis- 
tances with log-log scales. N 3 = 16K (k§ = 4.090 and k z = 2.200). d = 9, 10 and 
11. 




Figure 5: Genus distributions of mother boundary at the critical point for various dis- 
tances with log-log scales. N 3 = 32K («§ = 4.195 and k 3 = 2.222). d = 9, 10, 11 and 12. 
Solid line is drawn by the least squares of the data with d = 12 with 1 < g < 50 as a 
fitting range. 

The small mother boundary surfaces will suffer a large finite-size effect, so that almost 
all these surfaces become spherical surfaces. Therefore, we must introduce the lower limit 
iVcut) of the volume for the spherical mother surface in order to avoid the large finite-size 
effects. In this subsection we ignore spherical surfaces for simplicity. In next subsection 
3.2, we actually introduce the lower limit (V cu t) in order to measure coordination number 
distributions. Figs.|] and [5] shows genus distributions of the mother boundary surfaces 
P(g) with various distances and with N 3 = 1QK and 32K i. We naively expect that the 
genus distributions of the boundary surfaces will show the scaling properties and, in fact, 
we find scaling behavior , 

P(g) ~ g~ a , (5) 

where g is the genus and a is 0.55(1). a is obtained by the least squares from the data in 
Fig.|5] with N 3 = 32K and with 1 < g < 50 as a fitting range. It reveals that the scaling 
property of the genus distributions becomes the more clear the bigger the size of boundary 
becomes. Then, this scaling property will remains after the thermal limit N 3 — > 00. If Ko 
is out of the critical point this scaling relation disappears. 

3.2 Genus distributions in the strong coupling phase 

We also measure the genus distributions of the mother boundary surface in the strong 
coupling limit (i.e., k = 0) with various N 3 and distances. Fig|] shows the distributions 

§In Figs.|| and || we prepared about 400(^3 = 16K) and 100(^3 = 32K) configurations using the 
Monte Carlo method. We took 160 start points (simplices) for N3 — 16K and and 320 for N3 = 32K par 
a configuration. 
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Figure 6: Genus distributions of the mother boundary surface in the strong coupling phase 
for various sizes. 

of the mother boundary surfaces with three different volumes (N3 = 4K, 8K and 16K) 
0, and it reveals that these distributions seem to the Gaussian. In this measurement we 
select a distance (d m ) at which the peak value of each distribution becomes maximum, 
and actually in Fig.| we use the values d m = 9, 10 and 11 for N 3 = 4K,8K and 16K, 
respectively. The peak value (g p ) of each distribution becomes the larger the bigger the 
size of boundary surface becomes. We know that the manifold becomes crumple in this 
phase, and the fractal dimension seems to diverge. We, therefore, reasonably conclude that 
g p diverges when N 3 goes to infinity and that the manifold becomes crumple boundlessly 
in the strong coupling limit. 

3.3 Coordination number distributions of boundary surfaces 

We must look more carefully into these boundary mfds. In the last few years, several 
articles have been devoted to the study of boundary mfds. In two dimensions, it is 
revealed by ref. JHJ that the dynamics of the string world sheet (random surfaces) can 
be described by the time I evolution of boundary loops. Furthermore, the work |14| is 
also based on the idea that the functional integral for 3 and 4D quantum gravity can 
be represented as a superposition of less complicated theory of random surfaces 0. It is 
precisely on such grounds that we claim that the higher dimensional complicated theory 
of quantum gravity can be reduced to the lower dimensional quantum gravity. 

The question is that the closed and orientable two-dimensional boundary surfaces 
mentioned above can be recognized as the two-dimensional random surfaces described 

^In Fig.§ we have prepared 800(iV 3 = AK),400(N 3 = 8K) and 200(7V 3 = 16K) configurations and 
obtain 32, 000 genus data points evenly for each volume, 
"i.e., geodesic distances 

**In this case a definition of a time direction is different from our definition of a time slice, i.e., geodesic 
distance. 
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Figure 7: Normalized coordination number distributions for the mother boundary universe 
near to the critical point with A^ = 8K(k = 3.950 and k 3 = 2.165), 16K(k = 4.090 
and K3 = 2.200) and 32K(k = 4.195 and % = 2.222) with linear-log scales and with 
V&it = 200. We have measured the distributions only for the mother boundary surfaces 
with spherical (S 2 ) topology. Broken line is a guide line and denotes the asymptotic 
scaling relation: Pn 3 (<i) oc e~ ln s q . 

by the matrix model or the Liouville field theory. In this subsection, we focus on the 
distributions of the coordination numbers. The distributions of the coordination numbers 



have been calculated in ref. [15|. The probability distribution Pn(q) can be extracted 
from the Green functions of the (ft 3 planar theory with spherical (S 2 ) topology without 
self-energy and tadpole graph in dual lattice, 

P N ^(q) - 16(-) > e , (6) 

where N is a number of 2-simplexes and q is the coordination number. The probability 
distribution Pn{q) with small q region is strongly dependent on the local lattice structures. 
The value (In 3) guarantees the randomness of the lattice (see Appendix A). In the pre- 
ceding subsection we pointed out that spherical mother boundary surfaces whose volumes 
were greater than V cut had to be considered. Fig.[7] shows the (correct-normalized) coordi- 
nation number distributions of the spherical mother boundary surface for various volumes 
with V cut = 200. In Fig.^we also show the result obtained by the theoretical calculation 
eq.(|6]) (broken line). The data of the coordination number distributions of the spheri- 
cal mother boundary surfaces become consistent with the distributions Pn{q) oc e~ ln s q 
for large enough q. On the other hand, there are a few discrepancies between the two- 
dimensional theoretical curve and the small q (q <~ 25) distributions in our simulation 
sizes. We cannot say from only these data whether the mother boundary surfaces with 
S 2 topology are equivalent to the random surfaces implemented by the matrix model or 
not. Further discussion will be presented in the next section. 



4 Summary and Discussion 



We investigate the boundary surfaces of three-dimensional DT manifolds with S 3 topology. 
It have been naively expected that the genus distributions of boundary surfaces will 
show the scaling properties near the critical point and, in fact, we find scaling behavior, 
P(g) ~ g~ a , where a is about 0.5 in our simulation sizes. On these grounds we may 
conclude that the thermodynamic limit (N 3 — > oo) can be taken in this scaling region. 

At the strong coupling limit (kq = 0) the genus distribution of mother universe seems 
to Gaussian distribution, and it is non-universal in the sense that the peak value of genus 
g p diverges as N 3 — > oo. 

On the other hand, in the weak coupling region, all of the size of the boundary surfaces 
are very small, ~ 0(1), and then these genuses are strongly restricted to zero (i.e., sphere). 
As a result, there is no scaling behaviour of the genus distributions in the weak coupling 
phase. 

Our numerical results show that the coordination number distributions of the spherical 
mother boundary surfaces of three-dimensional DT manifolds are consistent with the the- 
oretical prediction in the large q region. When these boundary surfaces can be recognized 
as surfaces of the matrix model or the Liouville theory, three-dimensional DT manifolds 
can be reconstructed by the direct products of X (two-dimensional DT surfaces) and d 
(geodesic distance). In order to confirm the equivalence between the boundary surfaces 
in three-dimensions and random surfaces in two- dimensions the string susceptibility 7 str 
of the boundary surfaces or the loop-length-distributions (LLD) of the boundary surfaces 
must be measured and compared with theoretical predictions. We have obtained some 
preliminary results of the string susceptibility exponents of the boundary surfaces, and a 
more complete numerical analysis will be published elsewhere. 

We will be able to apply our numerical analysis to other physical process such as the 
topology changing of surfaces of (2 + l)-dimensional gravity. In this process we substitute 
the geodesic distance (d) for the time, and the two-dimensional surfaces are recognized as 
the boundary of three-dimensional Euclidean manifold. 

Furthermore, the boundary surfaces of three dimensions in general have other nontriv- 
ial structures such as Hopf 's link and the nontrivial knot in the knot theory. Up to now, 
it is difficult to extract the informations for the extrinsic geometry, i.e., how to embed 
the boundary surface into R 3 because we have only intrinsic geometries of DT manifold. 
To know a mechanism of the formation of nontrivial links or knots is one of challenging 
subjects in the statistical mechanics. 



Appendix A 

We can obtain the asymptotic coordination number distribution: eq.(|]) only assuming 
the Poisson distribution function as follows, 

P(q) ~ e- aq . (AX) 

By using this and the lower limit of q > 3 we can calculate an average coordination 
number < q > 

Eqe^ = d a^oo 

Ee- ff? da v ^ 3 ; \-er° v ' 

On the other hand, < q > is trivial in two dimensional DT mfd, 

<9> =^= 6 -i" 6 - (A3) 



where x 1S the Euler number. From eqs.(A.2) and (A. 3) we obtain a = In g |16 
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